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ABSTRACT 

0^ ' Second order perturbation theory predicts a specific dependence of the bispectrum, or three- 
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point correlation function in the Fourier transform domain, on the shape of the configuration 
of its three wave vector arguments, which can be taken as a signature of structure formed by 
gravitational instability. Comparing this known dependence on configuration shape with the 
weak shape dependence of the galaxy bispectrum has been suggested as an indication of bias in 
the galaxy distribution. However, to interpret results obtained from projected catalogs, we must 
£C) • first understand the effects of projection on this shape dependence. 

We present expressions for the projected power spectrum and bispectrum in both Cartesian 
and spherical geometries, and we examine the effects of projection on the predicted bispectrum 
^ • with particular attention to the dependence on configuration shape. Except for an overall numer- 

ical factor, for Cartesian projection with characteristic depth D* there is little effect on the shape 
dependence of the bispectrum for wavelengths small compared to D* or projected wavenumbers 
qD* 1. For angular projection, a scaling law is found for spherical harmonic index £ ^> 1, 
but there is always a mixing of scales over the range of the selection function. For large I it is 
sufficient to examine a small portion of the sky. 
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■ 1. Introduction 



Subject headings: large-scale structure of universe 



The behavior of cosmological density fluctuations in perturbation theory is becoming a mature, well- 
understood subject, allowing analytic calculations of the expected behavior of correlation functions to high 
order (Fry 1984; Goroff et al. 1986; Bernardeau 1992; Jain & Bertschinger 1994; Scoccimarro & Frieman 
1996a, b; Scoccimarro et al. 1998). Starting with Gaussian initial conditions at early times, the first nonvan- 
ishing contribution to the n-point correlation function requires knowledge of perturbations to order n— 1. 
To leading order irreducible point or windowed moments of density follow the hierarchical pattern, 

Cn = (S n ) c = S n Q- 1 , (1) 

first found in observations, where the full galaxy three-point correlation function £ can be written as a sum 
of products of two-point functions £ as (Groth & Peebles 1977) 

Cl23 = Ql23 (£l2£l3 + £12^23 + £l3£23 )■ (2) 

In observations, of both the three-point function in space and the bispectrum in the Fourier transform domain 
(Fry & Seldner 1982), the dimensionless amplitude Q depends only weakly on overall scale or on configuration 
shape. In perturbation theory, the reduced three-point function amplitude, or the corresponding bispectrum 
amplitude in Fourier space (see eq. Jl6| below), is not a constant but a function of the lengths of the triangle 
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or of the three wave vectors that make up the configuration (Fry 1984; Bouchet et al. 1995), a function that 
varies with the power spectrum index n but depends only weakly on cosmological parameters such as the 
fraction of critical density fi an( A the cosmological constant A (Bouchet et al. 1992; Bouchet et al. 1995; 
Scoccimarro et al. 1998). The apparent disagreement between theory and observation can be understood as 
an effect of bias. A nonlinear but local bias in the galaxy distribution, with fractional contrast in the galaxy 
number density determined as a function of the local mass density contrast, 6 g = J] bkS k jk\ rj 6(5 + ^&2<5 2 , 
gives a galaxy three-point amplitude that depends on the first two coefficients in the expansion (Fry & 
Gaztahaga 1993), 

Q 9 = \Qi2 3 + ^, (3) 

where b = b% is the usual linear bias factor. Matarrese, Verde, & Heavens (1997) have estimated the expected 
uncertainties and analyzed numerical simulations to investigate the practicality of applying the method to 
data, and have concluded that the method works well in the unbiased case, with likelihood centered around 
b = 1 and 62 = 0. The local bias model may be an oversimplification, but represents the minimal nonlinear 
degree of freedom in any nonlinear galaxy formation model. Comparison of the observed galaxy three-point 
amplitude in the Lick catalog with that expected in perturbation theory using equation (^ provides an 
estimate for b that is larger than expected from other methods (Fry 1994). 

Any measurements that provide information about bias are important for our ability to infer proper- 
ties of the cosmological mass distribution from galaxy observations. However, many of the largest galaxy 
compilations, including the Lick catalog, are angular positions only, and an important question is whether 
angular projection has a significant effect on the results. For angular correlation functions and moments this 
has been studied for some time (Limber 1953; Rubin 1954; Groth & Peebles 1977; Kaiser 1992; Polio 1997; 
Gaztanaga & Bernardeau 1998). In this paper we compute the effects of projection in the transform domain, 
on the power spectrum and bispectrum. In the following, in Section 2 we first compute a simple Cartesian 
projection, where we can easily observe the approach to the asymptotic small and large scale results, and 
present numerical results in Section 3. Following this, we examine angular projection in Section 4, over the 
full sky and also over a small patch. Section 5 contains a final discussion. 



2. Cartesian Projection 

For simplicity, we first present results for projection along a Cartesian direction, denoted z. Cartesian 
projection misses an important scale mixing found below in projection at fixed angular positions, but it is 
simpler both to compute and to understand, and thus provides a useful guide in understanding later results. 
The projected number density of galaxies n p is the space density n(x,y,z), summed over z, weighted by a 
selection function F(z) that gives the probability that an object at z is included, 

n p (x,y) = J dz F(z)n(x,y,z). (4) 

For convenience, although this is not necessary we take F to scale with a typical depth D* , F — F(z/D*) 
(Peebles 1980). In an attempt to separate projected and unprojected quantities, in what follows we will write 
projected functions with a subscript p: n p , £ p , etc. To distinguish positions in the two-dimensional projected 
space from those in the full three dimensions, we will denote the former as x and the latter as r; similarly, 
wave vectors in the two-dimensional space we will denote q and in three dimensions as fc. It will occasionally 
be necessary to refer to a three- vector made from a two- vector and a third component; r = (x, z), k = (q, fc), 
etc.; k — q means k = (q, 0). Thus, equation (|J) can equivalently be written n p (x) — J dz F(z) n{r). 
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We wish to study statistics of the projected density. In space, we write (n(r) ) = n, and the two-point 
correlation function £ and the three-point function £ are defined by 

(n(n)n(r 2 )) = n 2 [1 + £(r 12 )] (5) 

and 

( n(n)n(r 2 )n(r 3 ) ) = n 3 [1 + £(r 12 ) + £(r 23 ) + £(r 31 ) + £(r lf r 2 , r 3 )] , (6) 

where angle brackets ( • • • ) indicates an ensemble average, equivalent by homogeneity to an average over 
space, and ri 2 is shorthand for r\ — r 2 . Cosmological symmetries imply that £ and £ depend only on 
distances between points. 

Similar expressions hold for projected functions. The average projected density is n p = j dz F(z) n, 
while from equation (||) we can express moments of n p as integrals over the spatial correlation functions, 

n 2 p ti p {x 12 )=n 2 J dz 1 dz 2 F(z 1 )F(z 2 )Z(\r 12 \), (7) 

n 3 (p(xi,x 2 , x 3 ) = n 3 J dzidz 2 dz 3 F(zi)F(z 2 )F(z 3 ) £(ri, r 2 , r 3 ). (8) 

These are equivalent to relations for projected angular functions (Limber 1953; Peebles 1980). 

Because correlations depend only on separations, it is natural to change variables from zi, z 2 to the 
appropriate mean z = \{z\ + z 2 ) or z = \{z\ + z 2 + z 3 ) and differences z X2 — Z\ — z 2 . In the approximation 
that correlations are small except for small separations, we can write 

fdzF 2 (z) r 
f dz F 3 (z) f 

C P (x 12 ,x 13 ) = ^ dzF ^Y J dz 12 dz 13 ((r 12 ,r 13 ) (10) 

From these, it is apparent that if the spatial three-point function is hierarchical with constant Q (eq. [||), 
then so is the projected three-point function, with 

Q p _ [JdzF 3 (z)][JdzF(z)] 

Q [jdzF 2 {z)} 2 ■ [ 1 

This result is independent of the normalization of F. By the Schwarz inequality, Q p /Q > 1. The equivalent 
result was obtained in the angular case by Fry & Seldner (1982). 

Our goal now is to compute projected moments in the Fourier domain. In Fourier transform space we 
have the amplitude 

n(k) = J d 3 rn(r)e- lk r , (12) 

n(r) = J-^n(k)e ikr . (13) 

The mean density contributes (2Tr) 3 S-£>(k)n to h{k). Moments of the Fourier amplitude for k ^ give the 
power spectrum P(k) and the bispectrum B{k\ 1 k 2l k 3 ) (the Fourier transforms of £ and £ respectively), 

(n(fci)n(fca)) = n 2 [(2 7 r) 3 fe(fe 1 + fc 2 )] P(fc), (14) 
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(n(fci)n(fc 2 )fi(fc 3 ) > = n 3 [(2 7 r) 3 (5 D (fc 1 + k 2 + k 3 )} B(k u k 2 ,k 3 ). (15) 

The "momentum conserving" Dirac 5-function, arising from homogeneity, ensures that the bispectrum is 
defined only for configurations of wave vectors that form a closed triangle, ^2 = 0- The reduced three- 
point amplitude is defined to be 

n( , , , x g(fcl,fc 2 ,fc 3 ) 

g( 11 2 ' 3) ~ P(fc 1 )P(fc 2 ) + P(/fci)P(fc 3 ) + P(fc 2 )P(fc 3 ) (16) 
The projected Fourier amplitude is 

= | <frf (*) | ^Lft(*,)e<».-'e-'»" = j ^#*(*,)8(«,M. (17) 



the spatial amplitude evaluated at fc = q smeared in the projection direction over F(k z ) = J dz F(z)e lk * z . 
Note that if F has typical scale D*, then F has extent 1/D*. 

Moments of the projected amplitudes then give the corresponding projected spectra. Normalized by n p 
these are 

™ = ]j^Jtr^' p M +k ^ < 18 > 

Bp(q 1 ,q 2 ,q 3 ) = [j dzF ( z ^3 J ^T^T F ( k zi)F(k z2 )F(~k zl - k z2 ) B(q 1 , k zl ; q 2 , k z2l q 3 , k z3 ). (19) 
These are the basic results of this section. 

The large- and small-gP^* limits of these expressions are both of interest. Analogous to the small 
separation approximation in coordinate space, when D* is large, F is nonvanishing only for small k z ; the 
arguments of P and B are then the transverse q and can be taken outside the integrals. In this limit the 
projected spectra are, up to a constant factor, the corresponding spatial function evaluated in the transverse 
plane, k « (q, 0). By Parseval's theorem and a simple generalization, the integrals over k z can be written as 
the same integrals over z that appear in equations @j and (|To|), 

™- p ^{0W' ,20 > 

f dz F 3 (z) 

B p (qx,q 2 ,q 3 ) -*■ B( qi ,q 2 ,q 3 ) jj dzF ^ 3 , ( 21 ) 

x \ [fdzF 3 (z)][fdzF(z)} 

Q P ( qi ,q 2 ,q 3 )^Q( qi ,q 2 ,q 3 ) U { J dzF2 \ z)]2 ( 22 ) 

(qD* ^> 1). The opposite limit, large wavelength or small D* , corresponds to data confined to a thin sheet 
or slice taken out of three-dimensional space. In this limit, we have F(k z ) w F(0) = J dzF(z), and 

/ ^Pl(<? + k 2 z ) 1/2 ] (23) 

B p (q 1 ,q 2 ,q 3 ) -> J B(q x , k zi ; q 2 , k z2 ; q 3 , k z3 ) (24) 

(qD* -> 0), where k z3 = -(k zl + k z2 ). 
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3. Numerical Results 

In this section we presents results of numerical integrations of equations ( |i^ , |l9| ). For numerical calcula- 
tions we take F(z) = (zj 'D*)e~ z l D , F(k z ) = 1/(1 — ik z D*) 2 . This analytically convenient form mimics the 
relative lack of nearby galaxies because of geometry and the rapid fall of the luminosity function above L* . 
Because integration is a smoothing operation, we expect the general nature of the results to be insensitive 
to the precise behavior of F. For this F, integrals appearing above can be done analytically, 

I n = [ dzF n (z) = n! (D*/n) n+l , (25) 



(26) 



and so from equation ( |22[ ) we have for large qD* 

Qp _ Iila _ 32 
~Q ~ If ~ 27' 



For our numerical study we use the perturbation theory form for the shape dependence of B, 

B(k u k 2 , fc 3 ) = Q(k 1 ,k 2 )P(k 1 )P(k 2 ) + Q(ki,k 3 )P{ki)P(k 3 ) + Q(k 2 ,k 3 )P(k 2 )P(k 3 ) (27) 



10 , t t , (ki . kA . 4. 



Q(k i ,k j ) = f + (ki ■ kj) I f + '-f J + j(k t ■ kjf (28) 



where 



(Fry 1984). 

Figure [l] shows the projected P p (q) (eq. for n = +1, 0, — 1, and —2, evaluated for D* = 1. Results 
scale with D* ^ 1 as P p (q; D*) = D*~ (n+1) P p (qD*). Dashed lines show the large-g result (eq. j20|), P p - q n . 
There is a change of characteristic behavior near q = 1. The behavior at small q, corresponding to scales 
much larger than the depth of the projection, is that of data on a thin two-dimensional plane sheet or slice 
embedded in three dimensions. 

Figure ^ shows the dependence of the projected three-point amplitude Q p (0) on configuration shape for 
triangles with sides q\ — q, q 2 = q/2, separated by angle 9. Curves are shown for qD* = 0.001, 0.01, 0.1, 1, 
10, 100, and 1000, as indicated in the caption. The four windows show results for power spectrum indices 
n = +1, 0, —1, and —2. In some cases, the curves for small qD* overlap and are indistinguishable, as is also 
always the case for large qD* . In the l&rge-qD* limit the projected result is the result in space multiplied 
by an overall factor, here 32/27. 

Projection changes the shape of of Q, making the curves in Figure ||| flatter at shallow depth, or at 
large wavelength, just as would a bias with b > 1. Misconstruing the flattening as an effect of bias gives the 
parameters b = b\ and b 2 /b 2 listed in Table 1. The rms difference between the projected Q and the fit Q, 
A 2 = J d6 (Q p — Qb) 2 /^, is small, at most 0.0036, for n = —2 and qD* = 1 and often much smaller. The 
difference between the fit value of b and 1 is of order 1/qD* when qD* is large. 



4. Angular Projection 

4.1. Full Sky Coverage 

Typically, observations without distance information give galaxy positions not in Cartesian coordinates 
but in angle on the the sky; so we next consider the effects of projection onto the sphere. In spherical 



- 6 - 



coordinates, a position in space is specified by distance r from the origin and angular coordinates 6, <p, also 
denoted variously by solid angle and unit vector n. Summed over r, the projected angular density on the 
sphere is 

Af(9,<p) = J r 2 drF(r)n(r,6,<j)), (29) 

with average value N — n J r 2 dr F(r). 

The natural transform on the sky is by spherical harmonics. For comparison with later results, we 
compute first the full spherical harmonic transform. First, we compute the amplitude 

fiflm = [ dQM{n)Y e * m . (30) 



We wish to to relate moments of this amplitude to the spatial power spectrum. The first step is to express 
Afim in terms of the Fourier amplitude n(k) defined in equation ([l2|). Useful relations towards this end 
include the expansion of a plane wave in spherical coordinates, 



e ikz 



= JV(2* + l)j<(fcr)P*(cos0) (31) 



and the addition theorem for spherical harmonics, 

47r^F, m (n 1 )y/ m (n 2 ) = (2t+l)P t (h 1 • n 2 ), (32) 

m 

where Pg is the Legendre polynomial of order £ (Jackson 1975). From these, we obtain the transform 
amplitude 

/d 3 k 
—-^h(k)F t (k)A7ri e Y; m (k), (33) 

where Fg{k) = J r 2 dr F(r) je(kr) is real. 

Only TVoo has a nonvanishing expectation, (Moo) = V^nFo(0). The second moment ( Afe 1 m 1 Ne 2 m 2 ) 
vanishes unless i\ = £2, mi = —fn%. Normalized by ( A/qo ) we have 



\M^ m \ 2 ) = f d 3 k Ff{k) 

(M) ) 2 J (27T)3^( ) 



Pt = ™_/ = 1 P(fc)> (34) 



independent of m. 



When £ is large, the sequence of functions Fe(k) approaches a scaling function, Fi(k) ~ £ 3 ^ 2 F s (kD* /£), 
where F s (y) is peaked around y ~ 1 and has width Ay ~ 1, as shown in Figure ||. With a rapidly decreasing 
F(r), this can be understood as the main contribution to Fi(k) arising from the first peak of jt(x), which 
appears at kD* w £. With guidance from the small angle limit below, we find F s = {n/2) 1 / 2 y- 3 e- 1 / y , which 
agrees with the solid line so well as to be indistinguishable in the figure. Thus, for large £ we have the 
asymptotic expression 

v - pm ^Sw? f M) p( '- ,id ' ) - (35) 

For power law P(k) - k n , equation @ gives P s {£) - P(£/D*) - (£/D*) n . Results for Vt for power law 
P(k) are plotted in Figure [| the scaling result P s {£) is plotted as dashed lines. The behavior is similar to 
that seen in Figure [l], approaching the scaling result for large £. 
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The third moment of the amplitude in equation (33) is 



(27r) 3 (27r)3 (27r)3 FiMFiMFiM {Airnfi^^ 

x[(2^) 3 <5 D (fc 1 + fc 2 + fc 3 p 123 r; imi (fci)y; imj (fe2)y < ; m! (fe3). (36) 

Indices ^i, £2, and ^3 specify the overall frequency of angular oscillation of Ye m , while mi, m,2, and m3 specify 
the orientation of the oscillation with respect to a chosen polar axis. Because the universe is believed to be 
statistically isotropic, we are interested in the part of this expression that is invariant under rotations. This 
can be expressed using properties of the Wigner 3-j coefficients (see Edmonds 1957), related to Clebsch- 
Gordan coefficients by 

= 77T. , tm/ 2 \3i mij 2 m 2 \j 3 -m 3 ■ 37) 
mi m 2 m 3 J (2j 3 + I) 1 / 2 

Under a rotation D(us) with Euler angles u> = (a, /?, 7), states with the same j but different values of m mix, 

D{u)\jm) = Y, V ™m'\jrn') (38) 

where the T> are matrix elements of D(u), 

vlil> = (jm'\D(u>)\jm), (39) 

with 



= ( _ ir ( ) 7, m (/3,a). .40) 



4tt 



1/2 



In terms of the rotation matrix elements, the effect of a rotation on the 3-j coefficient is 

\' pO'i) r,(j2) pOa) f ii J2 J3 \ = / 31 h h \ . / 411 
S™i ™2™2 m' 3 m 3 I TO ' 772,' m > \ mi m, mW ' 1 j 



t j '"'2 ) 3 

if we take mi = m2 = TO3 = 0, we then have 

1/2 f£i £ 2 I, 



V ( * 1 t2 h \ Y Y Y 
^ \ mi m 2 m 3 i 



(24 + 1) (% + !) (24 + 1) 



3 





(42) 



(4tt) 3 

Nonvanishing contributions to the sum have mi + m 2 + mz = 0; the must also obey triangle inequalities, 
\ti — £ 2 \ < £3 < £1 + £ 2 , etc.; and the sum £\ + 1 2 + £3 must also be even. Removing the dependence on m 
leaves a result independent of orientation of the polar axis. This leads us to consider the averaged angular 
bispectrum, 



mi : m 2 ) m 3 



h h £3 \ {J^i imi JVe 2m2 JVe 3 m 3 ) 



mi m 2 m 3 J (7V o) 

/■ d 3 *! d 3 k 2 FtM F t M F t3 (k 3 ) 

= J-f,fof-, / rr ~; rrr — = = = tS(K\, K 2 , real, 

Ws i (2tt)3 (27T) 3 Fo(0 ) f (0) F (Q) k ^ 

where k 3 = — fej — &2 and the coefficient I^tata is 



(43) 



/, 1 ^3 = [(24 + l)(2^2 + l)(24 + l)] 1/2 ^ 1 ^ £ Q 3 J. (44) 

This angular averaging was introduced by Spergel & Goldberg (1998) for the microwave background. This 
is formally very similar to the Cartesian result, equation ([l9]). For the behavior at large £ we look at partial 
sky coverage in the small angle approximation in the next section. 
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4.2. Partial Sky Coverage 

In practice, there is no galaxy or cluster catalog with complete sky coverage, and in deep catalogs 
interesting scales cover only small angles. It is far more efficient computationally to treat a small square 
patch of sky as flat and use a Fourier transform, or especially FFT, rather than the spherical harmonic 
transform (Fry & Seldner 1982). Our final goal is to determine just what the result is of such a procedure. 

The Fourier transform in "flat" angular coordinates 6 is 

n p (K) = [ d 2 9n p (0)e- lKS . (45) 
Ja 

where 9 = {9\ ,9 2 ) = (0 cos cj>, 6 sin (f>) runs over a small patch of sky say a square of side A, and k is a Fourier 
frequency conjugate to 6. The projected amplitude n p (K) in terms of fi(k) is then in full 



«) = / dHe~^ e I ^F,(fc)n(fc)^4^/y £m (n)y/ m (fc), (46) 



For the full sky transform, the integrals over angles involved orthogonal functions, leading to the simple 
result in equation (|33|). The projected second moment of n p here is 

00 

{h p {Ki)n p (K 2 )) =n 2 p J2A e {K 1 ,K 2 )V e (47) 
where Ve is as in equation (|34|) and 

A<(Ki,Ka) = (2*+l) / d 2 e 1 d 2 e 2 e- lK ^e- lK2 ^P e (h 1 ■ h 2 ). (48) 

J A 

This integral behaves as an effective ^-function, peaked for Kx + K2 = and at i = \k\, but smeared in k and 
I over roughly ±(2n/A), or ±1 wave per box. This is shown in Figure || for k — (8, 0), (7, 4), and (5, 6), in 
units of 27r/A = 18.947, each with approximately the same magnitude, k = 8, 8.062, and 7.810, respectively, 
but different orientations. Both the location and width of the peaks behave as expected. 

Useful expressions can be obtained when A is small, so that r = (r9 cos <f>, rO sin <fi, r). Integrals over 
angles such as 

+A/2 d6 x e 1 ^-^ = sin [ A ( fc * r ~ *^)/ 2 ] M 9 ) 
-A/2 X (k x r-K x )/2 

are effectively <5-functions, with resolution in k of ±27r/A. With this, the projected amplitude becomes 



h p (K) = I — [ drF(r)e tkr h(K/r,k)= f—f ^ F(l/u) e iKX n(nu, kuX), (50) 
1 2tt I 2tt I u z 



where r = 1/u and k = Xku. Then we have for the second moment 

(^(kxKM) =n 2 (2^) 2 <5 D (0 1 -^)- / ^ e iAA " / F(l/ Ul )F(l/u 2 ) P[k iUi (1 + X 2 ) 1 ^ 2 ], (51) 
' K 2 J 2tt J u{u^ 

where the (5-function in the second moment of n(k) require that K\ and k 2 are in opposite directions and 
thus u 2 k 2 = u\K\. For power law P(k) ~ k n the power spectrum factors to (1 + X 2 ) n / 2 P{niui). The integral 
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over A can be expressed in terms of modified Besscl functions of order ±(n + l)/2 for non-integer n, while 
for negative integer n, 

{(5 d (Ak) n = Q, 

$K (\Ak\) n = -l, 

i e -l A «l n=-2, ' 1 j 

ilAKl^dA/cl) n = -3. 

All of these are appreciable only for |Are| <; 1, fall off exponentially for large An, and integrate to 1; so 
for any n the integral is once again effectively a (5-function of An, to an accuracy ±1, where for small A, 
K is large. Requiring k 2 = M implies u 2 = u\ also, and thus the second moment of the projected Fourier 
amplitude in equation (pll) becomes 



( h p (Ki)np(K 2 ) ) = n 2 [(2tt) 2 (5d(ki + k 2 )] P p (k), (53) 

where 

p ^ = wm^ I^ F2ll/u]p(KU) - (54) 

Although this expression was obtained for a power law -P(fc), to the extent that the integral over A in 
equation (52) is effectively a (5-function, this ought to apply for any well-behaved power spectrum. This 
equation is the same Kaiser's (1992) equation (A. 2), obtained here by a different path, and is equivalent to 
the scaling result in equation ( |35| ) for P s (£) at large £. For F(r) — exp(— r/D*), this gives us the form of the 
scaling function, F s (y) = (ir/2) 1 / 2 y~ 3 e~ 1 / y . 

For the third moment, we have 

i~t \~ t \~ i w _ 3 fdk 1 dk 2 dk 3 f du x du 2 du 3 , 

J Z>7T ZiTT ZiTT J U-^ 

x [(27r) 3 (5 D (fci +k 2 + fc 3 )] B(uKi,ki;uK 2 , k 2 ;un 3l fe 3 ), (55) 

where as before u = l/r. One component of the momentum-conserving (5-function constrains k\ -\-ki~\-kz = 0; 
the two others ^2 UiKi — 0. One of these last, in any direction not perpendicular to K\, can be taken to give 
the value of ui, and the second component to give u 2 , in terms of u 3 and the Ki, leaving 

J ZTT J ZTT J ufllljiU^ 

x F(l/u\)F(l/u 2 )F(l/u 3 ) B(uiKi,u 2 k 2 ,u 3 k 3 ) (56) 

As was found for the second moment, the fc-integrals are effectively 5-functions, enforcing constraints that 
require the triangle in the Ki to be closed; when these are satisfied all three Ui are equal. Thus, we have 
finally 

(np(Ki)np(K2)np(K 3 ) ) = n 3 [(27r) 2 <5 D (Ki + k 2 + k 3 )] B p (ki,k 2 ,k 3 ), (57) 
B p ( Ki ,k 2 ,k 3 ) = { j duF ( 1/u y ui]3 J -^F 3 (1/u)B(uk 1iU k 2 ,uk 3 ). (58) 

This is what we have been aiming for, an expression for the projected bispectrum in terms of the full 
bispectrum in space and the selection function F. For power law P(k) ~ k n , B(uKi) = u 2n B(Ki), and 
so the projected bispectrum has the same shape dependence as the bispectrum in space, with an overall 
multiplicative factor 

_ [f duF^l/u)!! 2 ™-^! duF{l/u)u-^ 



[J duF 2 (l/u) u" 



'-Qua. (59) 
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Written in terms of r — l/u, this is the same projection factor obtained by Groth & Peebles (1977) for Q in 
space and by Bernardeau (1995) and Gaztahaga & Bernardeau (1998) for the skewness S3. 

As is apparent in equations ( |54| ) and (p^), at fixed angles we are probing clustering on physical scales 
A?* = rA8; and as r ranges over a broad selection function this can be a substantial variation. For a power 
spectrum that is not a simple power law, the effect of projection depends on how much the effective index 
n e gt = d\ogP/d\ogk varies over the extent of the selection function. An example is shown in Figures ^| and 
@ for the cold dark matter (CDM) power spectrum, as parametrized by Bardeen et al. (1986), for Q = 0.3, 
h = 0.7. The CDM spectrum has n = +1 as k — > and n — > —3 as k — > 00, and the projected spectrum has 
these limits also. Where the effective index is almost constant, the projected power is the same as would be 
obtained for a power-law spectrum with the same n. However, for values of k near the peak the effective 
index varies enough to give a visibly different result. For the bispectrum, the effects of projection include 
a change in the amplitude of the variation with angle, the effect of a linear bias parameter 67^ 1; or as 
an overall shift in value, similar to a quadratic bias parameter 6 2 7^ 0; or a change in the dependence on 
configuration that is not well fit by a local bias. Over the range of scales of the Lick catalog (Fry 1994), 
which corresponds to n = 25-150, the apparent bias parameters in Table 2 differ from 1 by about 10% to 
only 2% at the scales most heavily weighted. 



5. Discussion 

In this paper we have investigated the effects of projection on the power spectrum and bispectrum in a 
variety of geometries. In general, projection is a smoothing operation, smearing out the dependence of the 
power spectrum and bispectrum on scale and the dependence of the bispectrum on configuration shape. 

The degree of smoothing depends on the depth of the projection. For projection along a Cartesian 
direction we observe a range of behavior as this depth varies from very shallow to very deep. On scales 
small compared to the depth or wavenumbers q such that qD* is large, except for an overall numerical factor 
projection has minimal effect; but the result of smoothing becomes increasingly visible as q decreases, with 
a qualitative change in behavior at qD* « 1. For power law P(k) ~ k n , the projected power spectrum is 
also P p (q) ~ q n with the same value of n for large qD* but the slope becomes noticeably shallower for small 
qD* . For the bispectrum we can quantify the result by the apparent bias factors given in Table 1. For large 
qD* the apparent bias parameters differ from 6=1 and 62 = by an amount of order 1/qD* or smaller. 



This can be traced to equations ( |18[ ) and (19), where it can be seen that the range of scales contributing has 
spread Ak/k ~ 1/qD*. 

The projected full sky angular distribution has power Vi in spherical harmonic amplitudes as given in 
equation (gjj). For a power-law power spectrum, the scaling limit in equation ( p5[ ) gives Vt — > P s {£) ~ t n 
with the same n, while the exact result at small I departs from this scaling limit in a way similar to the 
behavior of Cartesian result as qD* — > 1. Restricting the angular transform to a small portion of the sky of 
size A smears the result over values of I over a range ±27r/A. Once again the limiting behavior for large k, 
where k is a transform variable conjugate to 9, is P p (n) ~ K n . 

Additional effects enter when the power spectrum is not a pure power law, The integrals in equation ( |54| ) 
and equation ([58J), and the obvious generalization for the n-point spectrum 

1 f du 

B n , p (n 1: ...,K n ) = ^ duF{1 / u)/ui]n J ^F n (l/u)B n (uK 1 ,... )U K n ) (60) 
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Ki = 0), average the spatial function over the range of scales covered by the selection function F(r). When 
the effective index n e ff = dlogP/dlogk has little dependence on scale, the power-law results still remain 
approximately valid, but when n e s is changing appreciably over the scales sampled by F the projected 
functions can vary greatly, manifesting as apparent values b and b 2 - 

From what we have learned above, to avoid introducing artifacts of projection in analyses of data we 
should limit consideration to scales where k ^> 1, and for angular transform on a domain of side A to kA > 1 
as well. These restrictions appear to be satisfied in the Lick analysis (Fry 1994), for which k = 25-152 and 
kA = 31-188. This would appear to suggest that the weak shape dependence of the Lick bispectrum (Fry 
& Seldner 1982; Fry 1994) cannot be entirely ascribed to projection (but could be from nonlinear evolution, 
Scoccimarro et al. 1998). 

Planned future microwave background experiments will eventually lead to precise measurements of cos- 
mological parameters such as the Hubble constant Ho, the fraction of critical density Slo , and the cosmological 
constant A, and also of properties of the primordial power spectrum such as spectral index n and normaliza- 
tion Q. On the scales of galaxy clustering, what we will then be studying will be details of the modulation 
of,primordial power by the dark matter. In order to do that, we must be able to measure bias parameters. 
Study of the galaxy bispectrum appears to be one tool with the potential to do this, and projection does 
not appear to introduce unsurmountable difficulties. 

A portion of this work was done while J.N.F. was a visitor at the Aspen Center for Physics. A similar 
subject has been addressed recently by Buchalter, Kamionkowski, & Jaffc (1999). Research supported in 
part by NASA grant NAG5-2835 at the University of Florida. 
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Table 1. Cartesian Projection Bias Parameters 



Ti 


nil* 
qu 


1 lh 
1/0 


hn/h 2 


A 

u rms 


+1 


0.01 


-0.0015 


0.3598 


8.65 x 10~ 6 




0.1 


-0.0668 


0.4159 


4.18 x 10~ 4 




1 


0.2911 


0.4363 


5.37 x 10~ 4 




10 


0.9858 


0.0128 


1.14 x 10" 4 




100 


0.9993 


0.0004 


1.23 x 10" 5 




1000 


0.9997 


0.0003 


1.19 x 10" 5 





0.01 


-0.0023 


1.1689 


2.83 x 10~ 4 




0.1 


0.0031 


1.0335 


1.49 x 10~ 3 




1 


0.3605 


0.4826 


5.50 x 10~ 4 




10 


0.9493 


0.0358 


3.10 x 10~ 4 




100 


0.9984 


0.0013 


5.67 x 10" 5 




1000 


0.9993 


0.0009 


5.56 x 10" 5 


-l 


0.01 


0.0272 


1.4424 


9.04 x 10~ 4 




0.1 


0.0923 


1.0868 


1.48 x 10" 3 




1 


0.4519 


0.4341 


2.28 x 10^3 




10 


0.9511 


0.0361 


8.49 x 10~ 4 




100 


0.9984 


0.0017 


1.49 x 10~ 4 




1000 


0.9993 


0.0012 


1.46 x 10" 4 


-2 


0.01 


0.2097 


0.9742 


5.94 x 10" 4 




n i 

U. -L 


0.2521 


0.8458 


u7yi n~ 4 

u.l 1 A 1U 




1 


0.5250 


0.3788 


3.60 x lO- 3 




10 


0.9510 


0.0372 


1.33 x 10~ 3 




100 


0.9984 


0.0020 


2.34 x 10~ 4 




1000 


0.9993 


0.0016 


2.29 x 10~ 4 




Table 2. 


Angular Projection Bias Parameters 




K 




1/6 


b 2 /b 2 


^rms 


10 


-0.543 


1.081 


0.183 


1.34 x 10~ 2 


25 


-1.289 


1.019 


0.167 


8.17 x 10^3 


50 


-1.733 


0.993 


0.109 


3.84 x 10~ 3 


100 


-2.026 


0.981 


0.068 


1.37 x 10~ 3 


150 


-2.152 


0.978 


0.054 


6.62 x 10~ 4 


500 


-2.405 


0.977 


0.029 


3.53 x 10~ 4 


1000 


-2.499 


0.980 


0.022 


3.74 x 10~ 4 
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Fig. 1. — Projected P p (q) for n = +1, 0, —1, —2 evaluated for D* = 1. Dashed lines show the large-g limit. 
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Fig. 2. — Projected Q p (9) for configurations with q\ = 1, q 2 = 1/2, separated by 8, for power-law power 
spectra P ~ k n with n = +1, 0, —1, —2 as labelled. In each window, the solid curve shows the unprojected 
result scaled by the factor 32/27; the dotted curve shows the projected result for D* = 0.001; the dot-short 
dash curve D* = 0.01; dot-long dash curve D* — 0.1; short dash- long dash curve D* = 1; long dashed curve 
D* = 10. Results for D* = 100 and D* = 1000 are indistinguishable from solid curve. 
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Fig. 3. — Projected angular power Vt vs. £. Symbols show integrals over spherical Bessel functions for I = 
through 10. Filled triangles denote n = +1, open triangles n = 0, filled squares n = 0, and open squares 
n = — 2. Dashed lines show scaling result expected to hold for large t. 




Fig. 4.— Scaled £ 3 / 2 F e vs. k/l. Symbols show i = 10, dashed curve I = 100, solid curve I = 1000. The two 
curves are essentially indistinguishable. 
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Fig. 5. — Integral Ai(ki,K2) (eq. |Q) plotted vs. I. . Solid line shows K\ = (8,0) and k 2 = (—8,0), long- 
dashed line shows K\ — (7, 4) and k 2 = (—7, —4), and short-dashed line shows K\ — (6, 5) and k 2 = (—6, —5), 
in units of 27r/A with A = 19°. 
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1 



10 10 2 



Fig. 6. — Projected P p (k) obtained from equation (£4|) for the CDM power spectrum. The solid line shows 
the full integration of the CDM spectrum, while the long-dashed line shows what would be obtained for a 
power law spectrum with the effective local index n e s- Vertical lines mark the values k = 10, 25, 50, and 
150 used in the next figure. The dotted lines show the spread of scales that contribute to the integrand, 
k 3 x~ 3 F 2 (~k/x)P(x/D*) with x = u/kD*, for these four values of n. 
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Fig. 7. — Projected Q(9) obtained from equation (|5^) integrating over the CDM power spectrum for config- 
urations with K2 — 2 K i separated by angle 9 for four values of K = Ki, as indicated. The long-dashed line 
shows the actual integrated projection, the short-dashed line shows projection of a power-law spectrum with 
index n e ff evaluated at K\. The solid line shows Q before projection at scale k = k/D*, multiplied by the 
factor in equation (|59|). 



